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I*  IMHODOCTIOM 


Future  ballistic  aissiles  eay  need  the  ability  to  'spin 
during  boost  ascent  through  the  ateosphere  as  a  counter- 
■easure  against  first-generation  continuous-duty  laser 
threats.  since  the  conventional  missile  control  design  has 
uncoupled  pitch,  yaw  and  roll  channels,  some  performance  and 
even  system  stability  may  be  lost  in  rolling  the  missile. 
This  study  evaluates  the  performance  of  a  conventional 
pitch-yaw  controller  in  the  presence  of  spin.  A  conpensa- 
tion  scheme  is  then  proposed  to  restore  system  performance. 

The  block  diagram  of  a  submarine- launched  ballistic 
missile's  pitch- yaw  control  system  is  used  as  the  design 
vehicle  for  this  study.  since  the  physical  parameters  of 
this  system  are  classified,  only  the  block  diagram  is  used 
in  the  analysis  and  design  process.  All  numerical  values 
needed  for  analysis  and  simulation  are  borrowed  from  other 
systems  or  derived  as  needed. 

The  study  is  conducted  as  follows.  First,  the  basic 
laws  of  motion  and  thrust  vector  control  are  introduced. 
Together,  these  sets  of  eguations  define  the  flight  dynamics 
of  the  missile  and  the  method  for  controlling  the  missile. 
The  pitch-yaw  controller  for  the  non-rolling  missile  is  then 
designed  and  simulated.  This  system  has  uncoupled  channels 
for  pitch  and  yaw  control  which  become  coupled  when  the 
missile  is  rolled.  Next,  the  missile  is  forced  to  roll  and 
the  effect  cross- coup ling  has  on  the  system's  performance  is 
analyzed.  A  compensation  scheme  is  then  proposed  to  restore 
the  system's  performance.  Analysis  and  design  are  carried 
out  using  classical  and  modern  controls  analysis  and 
simulation  results  are  presented. 


1.  BISSILB  OIIABICS  ABO  IHBOSl  fBClOB  CCMIBOl 

Botation  and  Befecence  Axis  Systea 

Ibe  cefeience  axis  systea  used  tc  define  vector 
quantities,  diaensions  and  angles  of  interest  is  sboun  in 
Fig.  1-1.  Ibe  axis  systea  is  fixed  within  the  aissile  body 
with  its  origin  at  the  nissile*s  center  of  gravity.  The 
orthogonal  axes,  x,  y  and  z,  will  also  be  referred  to  as  the 
roll,  fitch  and  yaw  axes,  respectively.  Fig.  1-1  indicates 
that  the  right-hand  rule  is  used  to  deteraine  the  positive 
direction  for  rotational  vector  quantities. 


Fig.  1-1  Beference  Axis  Systea 


Table  1-1  lists  the  sjfsbols  used  to  cefec  to  kej 
parasetecs  and  quantities  of  interest.  Other  sjibols  sill 
be  defined  as  the;  arise.  The  aoaents  of  inertia  about  the 
center  of  gravity  (eg)  are  defined  as: 

A  =  2  tfs  (y**2*)  (1.1) 

B  =  Z  (z»»xa)  (1.2) 

C  =  2  (x»*y»)  (1.3) 

The  products  of  inertia  ace  defined  as: 

D  =  2  (^B)  J2  (1,4) 

E  =  2  («Jin)  *2  (1.5) 

F  =  2  {6'a)  xy  (1.6) 


Hhere  syBbolizes  the  sub  of  differential  aass  eleaents 

which  Bake  up  the  fissile  bedy. 


V.S'A*  s'A"  ■^■.  A."  ••'"  'V 


lABLf  1-1 

MOTATIOM  EEFINITIONS 


Roll  Axis 

Pitch  Axis 

' 

Yaw  Axis 

X 

y 

z 

Angle  of  incidence 

X 

0 

V 

Angular  rate 

p 

q 

r 

Component  of  missile 
velocity  along  axis 

u 

V 

w 

Component  of  force 
on  missile  along  axis 

X 

y 

z 

Moments  acting  on 
missile  about  axis 

L 

M 

N 

Moments  of  inertia 
about  each  axis 

A 

B 

C 

Products  of  inertia 

D 

E 

F 

- - 

2-  Eulei; *s  Equations  of  notion 

Ibe  equations  of  notion  foe  a  bod;  uitb  six  degrees 
of  freedon  and  a  lass,  n,  ace  suanacized  belon.  A  dot  (•) 
orec  any  sjobol  lefecs  to  the  tine  rate  of  change,  d/dt,  of 
the  quantity  indicated.  [Bef.  1] 


a(U  4  qw  -  rv)  =  X 


i(»  ♦  rU  -  pn)  =  X 


(1.7) 


(1.8) 


L«  »  •  ■ 


*  •'*  • 


n(M  -  qU  ♦  fv)  =  Z 


(1.9) 


Ap  -  {B-C)qr  ♦  C(r*-ga)  -  E(pq*r)  ♦  f  (rp-g)  *  1  (1.10) 


-  {C-A)rp  ♦  E(p*-r*)  -  F{9C*p)  ♦  D  (pq-r)  * 


-  {A-B)pq  ♦  F(qa-pa)  -  D|rp*q)  ♦  E(gr-p)  = 


(1.11) 


(1.12) 


Ihe  first  three  equations  above  are  force  equations 
and  the  second  three  are  noaent  equations.  In  particular, 
Eqs.  (1.11)  and  (1.12)  are  those  froa  Mhich  pitch  and  yaw 
angular  accelerations  are  found.  Mote  that  the  second  tera 
in  each  of  these  equations  constitutes  a  cross-coupling  tera 
between  the  two  channels  which  is  proportional  to  the  roll 
rate.  Ihis  cross-coupling  effect  is  eliaioated  in  the 
conventional  pitch-yaw  controller  design  by  preventing  roll 
and  thus  setting  the  roll  rate,  p,  equal  to  zero. 

Assuaing  that  the  aissile  body  is  approziaately 
rigid  leads  to  the  conclusion  that  the  products  of  inertia 
are  equal  to  zero,  thus  Eqs.  (1.11)  and  (1.12)  are  reduced 


Bq  -  (C-A)rp  =  B 


(1.13) 


Cr  -  (l-fl) pq  =  N 


(1.14) 


Thciist  Vector  Control 

Ihe  only  leans  used  to  control  the  flight  of  the 
■issile  under  study  is  to  deflect  the  exhaust  nozzle  as 
shown  in  two  dinensions  in  fig.  1-2.  reflecting  the  exhaust 
nozzle  causes  the  driving  force  of  the  rocket  notor  to  he 
pointed  in  a  direction  other  than  along  the  lissile's  longi¬ 
tudinal  axis.  Ihis  action  creates  force  conponents 
perpendicular  to  the  roll  axis  which  do  not  act  through  the 
■issile*s  center  of  gravity.  The  resulting  nonents  are 
about  the  y-  and  z-axes  and  cause  the  aissile  to  pitch  and 
yaw.  The  nonents  about  the  center  of  gravity  caused  by  the 


deflected  thrust  vector  are: 


H  =  (Tzllc 


(1.15) 


N  =  <V)lc 


(1.16) 


uhere  1^  is  the  length  fron  the  nissile*s  center  of  gravity 
to  the  point  of  application  of  the  thrust  (exhaust  nozzle) 
and  ly  and  T^  are  the  lateral  thrust  conponents  parallel  to 
the  y-axis  and  z-axis.  If  dg  and  d*^  are  the  angles  of 
nozzle  deflection  in  the  pitch  and  yaw  planes,  then  the 
perpendicular  thrust  conponents  are  related  to  the  total 
rocket  notor  thrust  by  the  following  eguations: 


ly  =  l'Sin(d^^) 


(1.11) 


Tj  =  T*sin(d0) 


(1.18) 


Fig.  1-2  Thrust  Vector  Control 


Assuniag  that  the  maxiaun  angle  of  nozzle  deflection 
is  snail,  as  it  oust  be  for  any  physical  systen  of  this 
iyp^r  pernits  further  sinplif ication  of  Egs.  (1.15)  and 


(1.16)  to  the  linear  Egs.  (1.21)  and  (1.22). 
siD(d0)  =  dg 

sin  (diu)  =  dy 


(1. 19) 


(1.20) 


n-21) 


N  =  I-dy-lp  n.'22) 

Equations  (1.21)  and  (1.22)  give  the  nonents  in 
pitch  and  yaw  caused  by  the  cocket  lotoc.  Coibining  these 
two  equations  with  the  general  rotational  nation  equations 
in  pitch  and  yaw,  Egs.  (1.13)  and  (1.14),  peraits  foraula- 
tion  of  the  general  differential  equations  for  pitch  and 
yaw. 

T-d0-l^  *  Bq  -  (C-A)rp  (1-23) 


T«d4;«lc  *  Ct  -  (A-B)pq  (1.24) 

Solving  for  the  angular  acceleration  conponents  gives: 
q  «  (1/B)(I‘d0-lg  ♦  (C-A)tp]  (1.25) 


r  =  (1/C)[l*dylc  ♦  (A-E)pq]  (1.26) 

Equations  (1.25)  and  (1.26)  ;;oint  out  the  cross'- 
coupling  effect  that  the  coll  rate,  p,  has  on  the  systen. 
For  instance,  if  the  roll  rate  is  set  equal  to  zero  the 
equations  for  pitch  and  yaw  ace  conpletely  uncoupled. 

q  »  (1/E)I-d0*l^  (1.27) 


c  *  (1/C)  T*  dy  Ig 


(1.28) 


Figure  1-3  deficts  the  rolling  pitcb-yav  sjsten  just 
described  in  blocJc  diagran  focn  and  highlights  the  ctoss- 
coapling  terns  introduced  b]f  roll*  Ihe  classical  nissile 
control  sjsten  design  approach  assuned  that  the  coll  rate 
vas  snail  and  uent  to  great  lengths  to  ensure  this  nas  true, 
flaking  the  assunption  that  the  coll  cate  is  negligible 
allows  the  system  designer  to  separate  the  pitch  and  jaw 
channels  and  deal  with  then  independentlj.  This  studj  will 
analyze  the  effect  cross-coupling  has  on  stability  and 
perfocnance  and  fornulate  a  conpensatlon  schene  for  the 
cross-coupled  systei. 


Fig.  1-3  Ibrust  Vector  Control  Dynanics 
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II.  DgggLQgHglT  Qt  tit  IQEr&OLLIIg  MISSILE 


The  strategic  eissile  whose  pitch-yaw  control  systea 
will  be  used  as  the  design  wehicle  for  this  study.  The 
block  diagraa  of  the  systea's  pitch-yaw  control  systea 
(without  roll)  is  shown  in  Fig.  2-1.  Since  the  physical 
paraaeters  of  this  aissile  are  classified*  only  the  block 
diagraa  will  be  used.  All  systea  constants  and  physical 
paraaeter  values  which  are  used  for  analysis  and  siaulation 
will  be  taken  froa  other  systeas  or  derived  as  needed. 

The  aissile  is  a  three-stage  booster  rocket.  Flight  of 
the  aissile  is  controlled  by  vectoring  the  exhaust  thrust  as 
described  in  the  preceding  section.  Based  upon  predeter- 
ained  flight  trajectories  and  aeasured  or  estiaated  paraae¬ 
ters  during  flight*  the  guidance  coaputer  generates  pitch 
and  yaw  coaaands  and  transaits  these  to  the  pitch-yaw 
controller  [lef.  2].  The  algoritha  which  the  guidance 
coaputer  uses  to  generate  the  pitch  and  yaw  coaaands  will 
not  be  addressed.  Only  the  response  of  the  pitch-yaw 
controller  to  a  general  set  of  coaaands  will  be  considered. 

Table  2-1  lists  the  physical  paraaeter  values  used  at 
the  various  stages  for  design  and  siaulation.  These  paraac- 
ters  bear  no  reseablance  to  those  of  the  Trident  II  aissile 
but  serve  only  to  provide  actual  nuabers  for  the  variables 
so  that  siaulation  could  be  done.  The  paraaeters  listed  in 
Table  2-1  were  actually  taken  froa  the  Saturn  V  Hocket;  also 
a  three-stage  booster  rocket.  The  aoaents  of  inertia  listed 
in  Table  2-1  were  obtained  froa  the  other  paraaeters  by 
assuaing  that  each  stage  has  a pproxiaa tely  the  saae  inertial 
characteristics  as  a  unifora  right-circular  cylinder. 


+  1 
re 


®c  + 


Nozzle 

Servo-Actuator 


Nozzle 

Servo -Actuator 


Fiq.  2-1  Noa-bolling  Pitch-Taw  Control  Systei 


Figure  2-1  shows  the  scheaatic  block  diagraa  of  the 
systea  used  to  control  pitch  and  yaw  in  the  aissile.  Since 
the  aissile  does  not  spin,  the  pitch  and  yaw  channels  are 
uncoupled.  The  systea  consists  of  a  servo-actuator  used  to 
deflect  the  exhaust  nozzle,  the  aissile  dynaaics  and  sensors 
which  feed  back  angular  velocity  and  position.  The  deflec¬ 
table  exhaust  nozzle  saturates  at  a  aaxiaua  aef lection  of  ♦ 
five  degrees.  In  this  analog  systea  all  gains  are  constant 
throughout  operation. 


TABLE  2-1 


FHESICAL  PABAHETEii  SPECIFIC AlIOM S 


1 

Blitial  Mass 

slugs 

1.488x10^ 

3.22x10^ 

8.14x10^ 

Mass  Hate 

slug/sec 

950 

81.49 

14.75 

Thrust 

lb, 

7.496x10^ 

1.125x10^ 

2.249x10^ 

Length 

ft 

135 

150 

100 

i 

■  . . i 

Burn  Time  sec  I5O  559  479 


A.  DESIGN  OP  TBE  SBBfO-ACIOATOB 


The  transter  functioa  of  the  3ecvo-dC  taatoc  which 
deflects  the  exhaust  nozzle  of  the  aissile  is  given  in  Lg. 
(2.1).  The  task  is  to  select  nueecical  values  foe  the 
pacaeeterSf  a,  b  and  f,  such  that  the  servo-actuator  has 
acceptable  perforsance  characteristics.  Again,  the  paraa- 
eter  values  will  be  selected  without  regard  to  those  used  on 
the  actual  eissile. 


d>  (s)  ^  _ abi _  (2.1) 

(s)  (s  ♦  a)  (s*  ♦  2fbs  ♦  b^) 

It  is  desired  that  the  actuator's  perforeance  reseable 
that  cf  a  well  daaped  second-order  sfstea  with  liniaua  over¬ 
shoot  and  a  short  settling  tine  in  response  to  step  inputs. 
These  perforaance  characteristics  need  to  be  aore  stringent 
than  those  iaposed  on  the  entire  systea  so  that  the  actuitor 
does  not  becoae  a  liaiting  eleaent  when  the  entire  systea  is 
considered  later.  The  design  specifications  used  for  the 
actuator  are  suaaarized  as  follows. 

Locations 

In  order  to  ensure  that  this  third-order  system 
behaves  essentially  as  a  second-order  systea,  the  real  pole 
must  have  a  magnitude  significantly  greater  than  the  real 
part  of  the  coaplex  roots  of  the  second-order  tera,  e.g. 

a  >  lOfb  (2.2) 

4 

2  .  Da  aping 

In  order  to  aake  the  doainant  coaplex  pair  of  poles 
exhibit  good  daaping,  the  value  of  f  aust  approach  one. 


3.  Settling  liac 

In  order  to  teep  the  actuator  troi  liaiting  the 
settling  tiae  o£  the  entire  systei,  the  servo's  settling 
tine  Bust  be  significantly  shorter  than  that  of  the  pitch- 
yaw  controller.  If  the  desired  settling  tiBe  of  the  entire 
systeB  is  to  be  approzinately  ten  seconds  and  the  actuator's 
settling  tiae  is  3X  of  the  systes's,  then  the  actuator's 
settling  tiae  will  be  selected  according  to  the  relationship 


below. 


’settling  servo 


<  0.03  t 


settling  system 


(2.4) 


4.  teak  Cwershoot 


Ihe  peak  overshoot  of  the  servo-actuator  aust  be 
snail.  Otherwise,  it  nay  have  an  appreciable  effect  on  the 
overall  systea's  peak  overshoot. 


n(servo)  -  l  < 


(2.5) 


Vhere  h  (servo)  is  the  peak  overshoot  of  the  actuator. 

Equation  (2.5)  peraits  selection  of  the  appropriate 
danping  ratio,  £.  If  the  systea  is  assuaed  to  exhibit 
second-order  characteristics,  then  Eg.  (2.6)  can  be  solved 
for  the  daaping  ratio. 


n  (servo)  -  1  =  exp 


’-  Iff  “1  <  0.01 


(2.6) 


Solving  for  the  daaping  ratio,  f,  yields  the  following. 


f  >  /  [ln(0.01)]g 
-v4*  ♦  [lii(O.Ol)  ]* 


=  0.826 


So,  a  daaping  ratio  of  0.9  will  be  used, 
also  satisfies  Eg.  (2.3). 


(2.7) 


Ihis  selection 


The  specification  of  tq.  12.4),  conbined  with  the 
above  choice  of  daapin9  ratio,  can  be  nsed  to  select  the 
natural  frequency  of  the  second-order  ten. 


^  settling  servo  “  j  ^ 


3  sec 


(2.6) 


Solving  for  the  natural  frequency  leads  to  a  selection  of  b 


lual  to  15. 


b  >  14.8  rad/sec 


(2.9) 


Equation  (2.2),  conbined  with  the  above  choices  of 
daaping  ratio  and  natural  frequency,  leads  to  a  selection  of 
the  real  pole,  a,  equal  to  ISO. 


a  >  lOffc  =  135 


(2.10) 


The  root  locus  of  the  systen's  characteristic  equa¬ 
tion,  obtained  by  varying  the  real  pole,  a,  is  shown  in  Fig. 
2-2.  Highlighted  ace  the  root  locations  for  the  above 
choice  of  paranetec  values.  As  shown  in  Fig.  2-2,  all  the 
coots  ace  well  into  the  left  half  of  the  s-plane,  a  large 
phase  naegin  exists  for  the  second-order  cowplez  pole  pair 
and  the  second-order  pair  is  far  enough  to  the  right  cf  the 
single  real  pole  to  achieve  doninance.  Figure  2-3  shows  the 
open-loop  frequency  response  of  the  actuator  and  aakes  the 
large  phase  nargin  aore  evident.  Figure  2-4,  the  closed- 
loop  frequency  response,  points  out  the  liiitation  of  the 
actuator  to  follow  input  signals  of  frequency  greater  than 
about  10  radians/second.  Figure  2-5  shows  the  step  response 
of  the  systes  and  as  indicated,  it  aeets  the  design  specifi¬ 
cations.  Note  that  although  the  sinulatioo  shown  in  Fig. 
2-5  is  that  cf  a  full  step  cf  unit  aagnitude,  the  nczzle 
actuator  saturates  at  a  laxinun  angle  of  deflection  of 
0.0873  radians  (5  degrees). 
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Fig.  2-5  Servo-Actuator  Boot  locus  Plot 


Figure  2-5  shows  that  the  aagnitude  of  the  real  pole 
is  far  greater  than  the  real  part  of  the  coaplex  pair  and, 
therefore,  will  have  little  effect  on  the  systei*s  perforw- 
ance.  Also,  the  doiinant  seccnd-order  pair  is  well  daiped 
with  approxiaatel;  a  72-degree  phase  nargin  and  an  effective 
daaping  ratio  of  0.95,  ihe  sjstea  is  stable  as  long  as  the 
tera,  ah*,  reaains  less  than  7.57  x  10». 
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Pig.  2-3  Servo- Actuator  Ofen-Loop  Frequency  Response 


21 


O’l  o»i-  o'gi- 


08C-  OI«-  ©’I 

(BP)  aaiuiNovn 


ou-  ••0«- 


Fig.  2-4  Ser»o-Actuator  Closed-Loop  Fceguency  Besponse 
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B.  0ES16B  OF  TBB  PITCB-IIB  COllBOLIBB 

Ibe  block  diagraa  of  the  pitch-]faH  control  systea  is 
shown  in  Fig.  2-1.  Since  the  pitch  and  yaw  channels  are 
identical  and  decoupled  only  one  of  thea  need  be  considered. 
The  paraweters  )  and  (1^4.)  will  fce  selected  so 
that  the  systew  will  eihibit  ainiaal  overshoot  to  step 
inputs  and  have  a  settling  tiae  less  than  ten  seconds.  The 
one-channel  systea  to  be  considered  is  shown  in  Fig.  2-6. 


-h  1 


(8+a)(8^+2fb+b^) 


la  ordec  to  tegin  the  analysis,  it  is  necessary  to 
tespocarily  ignore  the  nonlinear  saturation  of  the  servo- 
actuator.  Ihe  lineacized  transfer  function  fcr  the  pitch- 
yaw  control  systei  then  becones  that  of  £9.  (2.11). 


0(s)  =  KgabaTl^./E 

00  (s)  ss  ♦  s*(a*2£b}  ♦  s5|ba>2fbaj  ♦  s*  (ah*) 

♦  slK^gK^atail^/E)  ♦  K^ah*Hj,/E 


(2.11) 


Ihe  denominator  of  Eg.  (2.11),  when  set  egual  to  zero, 
defines  the  system's  characteristic  equation.  Ihe  charac¬ 
teristic  equation  nay  be  partitioned  as  follows.  Define  a 
new  variable,  J. 


J  =  at*R(Tl  /E) 
a  c 


(2.12) 


Ihe  characteristic  equation  now  becomes; 

ss  ♦  s*(a*2fb)  ♦  s3(b2*2ffca)  (2.13) 

♦  s*  (ab*)  ♦  s  (1^9  J)  +3  =  0 

Bearranging  yields  the  partitioned  form  shown  below. 


K^.s  ♦  1 


s*  *■  s*(a*2fb)  ♦  s3(b2  +  2fta)  ♦  s*(at*) 


=  -1 


f2. 141 


Equation  (2.14)  is  in  the  fora  for  root  locus  analysis 
where  J  is  treated  as  the  open-loop  gain  and  the  left  side 
of  the  equality  as  the  open-loop  transfer  function.  Varying 
Rre  and  J  together  sweeps  out  a  region  in  the  s-plane  into 
which  the  roots  of  the  characteristic  equation  can  be 
located.  Ihese  "root  relocation  zones"  are  shown  in  £ig. 
2-7.  [Bef.  3] 

Figure  2-7  was  made  by  varying  0  and  R re  in  Eg.  (2.14) 
and  plotting  the  root  locations.  tor  example,  if  Kre  =  0 
then  Eg.  (2.14)  contains  five  finite  poles  and  five  zeros  at 
infinity.  The  root  locus  will  follow  the  path  indicated  in 


Ay 


fig,  2-7  (as  J  increases).  note  that  the  real  pole  at  -150 
is  not  shown  in  Fig.  2-7.  As  Kr*  is  increased  incresentalljr 
fron  zero  towards  infinity  and  the  saie  procedure  is  used, 
i-e-,  plot  the  root  locus  over  a  large  (positive)  range  of 
J,  then  a  fasily  of  root  loci  will  be  swept  out.  Ihe  cross- 
hatched  area  in  Fig-  2-7  is  the  area  enlosed  by  the  root 
loci  for  •}  and  varying  fros  0  to  oo. 
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Figure  2-7  iodicates  that  the  syatea  is  unstable, 
two  coats  in  the  right  half-plane,  for  values  of  less 

than  or  equal  to  zero.  The  values  of  K and  J  ace  esti- 
■ated  as  0.8  and  1.57  1  IQs,  respectively.  This  selection 
locates  the  characteristic  egoation  coots  as  sbovn  in  fig. 
2-8.  The  open-loop  gain,  vas  found,  froa  the  definition 
of  J,  to  range  betveen  the  values  shown  below. 

1st  Stage  2nd  Stage  3rd  Stage 

1.071-0.07  1.51-0.  14  1.45-0.19 

Eiact  solutions  for  and  K were  not  obtained  analyt¬ 
ically  since  the  systea  contains  a  nonlinear  saturation 
effect  which  was  ignored  in  constructing  the  root  locus. 
Also,  K,  is  a  function  of  systen  paraaeters  which  change 
continuously  during  flight  and  dcanatically  at  booster  stage 
separations.  To  fine  tune  the  selection  of  and  , 

slnulation  studies  were  conducted  using  DSL/360  during  each 
of  the  three  booster  stages  of  operation.  The  results  of 
the  sinulations,  for  various  values  of  Rj.^  f  ate  shown  in 
Figs.  2-9,  2-10  and  2-11,  where  0  is  given  as  a  function 
tine. 

The  sinulations  illustrated  that  the  gain,  ,  has 

little  effect  on  the  systen»s  pecfornance.  This  is  because 
it  is  "upstreaa"  iron  the  saturated  nozzle  actuator  and,  as 
such,  can  only  shorten  the  tine  required  for  the  actuator  to 
reach  aazinum  deflection.  Since  this  tine  is  already  snail 
conpared  to  the  systen's  response,  has  negligible  effect. 
To  sinplify  the  probleo,  K,  was  set  equal  to  one.  The 
velocity  feedback  gain  has  a  pronounced  effect  on  both 

the  overshoot  and  the  settling  tine  of  the  systen.  Frcn  the 
sinulations,  was  selected  to  be  1.5.  This  choice 

ensures  the  specified  settling  tine  of  ten  seconds  during 
all  stages  and  results  in  overshoot  only  during  stages  one 
and  two,  where  it  is  less  than  ten  percent. 


/  f. 
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III.  IffBCl  OF  BOLL 


Rhen  the  noo-rolling  pitch-jaM  control  sjstea  just 
designed  is  forced  to  roll,  the  systen  fcecoaes  ccoss-coapied 
as  shown  in  Fig.  j-l.  The  lechanisi  which  creates  and  lain- 
tains  the  colling  lotion  will  not  be  considered.  It  is 
assumed  that  a  control  sjstei  outside  the  pitch-yaw 
controller  and  independent  of  the  lain  thrust  lotoc  lain- 
tains  an  appcoxiaately  constant  roll  rate,  p.  Furtherioce, 
it  is  assuied  that  the  roll  cate  can  be  Measured  during 
flight  or  estinated  icon  leasured  paraieters  with  reasonable 
accuracy.  The  purpose  of  this  section  is  to  analyze  the 
effect  roll  has  on  the  previously  non-rolling  pitch-yaw 
controller. 

The  Laplace  tcansforiations  of  the  linearized  non- 
rolling  systei*s  governing  differential  eguations  can  be 
written  as  two  separate  eguations;  one  for  the  pitch  channel 
and  one  for  the  yaw  channel. 

3*9  ♦  GK  j.gse  *00  =  G9g  (3.  1) 


s^r  ♦  HKy^st  ♦  HT  =  Hfg 
Nhere  G  and  H  ace  defined  as: 


(3.2) 


G 


_ ab* _ 

(s*a)  (s*  +  2fts*t*)- 


(3.3) 


H  =  K^aic/C) 


_ a^a _ 

-  (s*a)  {sa*2f bs*ta)- 


(3. ‘I) 
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Equations  (3.1)  and  (3.2)  can  be  rearranged  and  put  into 
■atriz  fora  as  shoen  in  Eg.  (3.5).  Mote  that  the  equation 
has  only  diagonal  non>zeco  eleaents  and  is,  therefore, 
uncoupled. 

s2«GKreS*G  0  9  _  G  <3]  [  ^  (3*5) 

0  s**HKr+s*HJ  uJ  [o  hJ  L'<c. 

a 

Rolling  the  zissile  leads  to  the  addition  of  cross¬ 
coupling  teras  between  the  equations  governing  pitch  and  yaw 
response.  The  transforaed  differential  equations  governing 
the  responses  of  the  rolling  systea  are: 

3*9  ♦  GKpg  s6  ♦  G0  ♦  QsV  »  GOg  (3.6) 


s*t  ♦  HKrysV  ♦  HV  ♦  BsS  =*  Etc  (3.7) 

tihere  G  and  H  ace  defined  in  Egs.  (3.3)  and  (3.4)  and  C  and 
B  ace  defined  in  Egs.  (3.8)  and  (3.9). 

C  =  (A  -  C)p/E  (3.8) 


F  =  (E  -  A)f/C  (3.9) 

Rearranging  Egs.  (3.6)  and  (3.7)  into  aatrix  fora  as  before 
leads  to  the  coupled  fora: 

sa*GKpg  3*G  ;s  ]9  _  [g  Olfflc]  (3.10) 

fis  s**HF^^s»hJW  ■  [o  rJ  [yc. 

Equation  (3.10)  shows  that  the  introduction  of  roll  has 
affected  the  transfer  functions  which  detecnine  the  input- 
to-output  relationships  of  the  systea  in  two  ways.  first. 


it  has  changed  the  direct  path  transfer  fnnctions  of 
pitch-to-pitch  and  law-to-yaw  connand-to^output.  Second, 
roll  has  introduced  pitch-to-yau  and  jaN-to-pitch  coniand- 
to-output  transfer  functions  ufaich  nere  rero  before  the 
■issile  vas  forced  to  roll.  Ihe  effect  that  these  changes 
have  on  the  systeB*s  stability  and  perforiance  uill  be  exan- 
ined  next. 

A.  EFFICT  OF  BOIl  01  SISIEB  SllBllITI 

As  pointed  out,  the  rolling  systea  is  governed  by  four 
transfer  functions  (2  inputs  x  2  outputs).  Consider  cnly 
the  pitch-to-pitch  input-to-output  transfer  function,  i.e,, 
hold  the  yaw-coaaand  input  equal  to  zero  and  ignore  the  yau 
output.  Under  these  conditions,  the  block  diagraa  of  the 
systea  nay  be  rearranged  as  shovo  in  Fig.  3-2.  Ihe  transfer 
function  obtained  froa  the  diagraa  is  given  in  Eg.  (3.11). 
As  indicated,  the  pitch  response  of  the  systea  is  influenced 
by  a  variety  of  systea  paraacters,  as  it  aas  before  the 
aissile  aas  rolled.  Nov,  however,  the  root  locations  of  the 
systea  are  also  affected  by  the  roll  rate. 


r.  w . 


9(S)  = 
0c  IS) 

where ; 


N^ss+N^S^^RjSa+NgS^^N^S^NQ 
s»  ‘•♦CgS*  ♦D3S">DyS»*DgS»*DjS**D^S* 
♦  Dj  sa^E^  S**C^  S*Cq 


(3.11) 


N5  *  (abVB)'IlcKe 

*  (aba/B)  31gK,(a42£b) 

Nj  =  (ab»/B)TlgKg(ha  +  2£ba) 

*  (abVB)IlcK«  (at*) 

Nl  »  (ab»/B)TlcF,(Rryl'yat2Tl^,/C) 
Nq  »  (afca/B)Tl^K^  (K^ab*Tlp/C) 

*  2(a*2fb) 
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Dg*  2(b**2£ba)  ♦  [  (A-B)  (i-CJ  p»/CB  ]  ♦  (a*2£t)2  =  Dg^  *  ^82^* 


Dy*  2ab*  ♦  2(a*2£b)  (b**2fba*[  <I-B)  {A-C)p*/CB])  =  D71  ♦  C72  p* 

Dg*  (b»»2fba)  (ba*2fba  ♦  2[  (A-E)  (A-C)  p*/CB ]) 

♦  (abaTlc)[  ]  ♦  (a*2ab)  (2ab*) 

♦  [  (A-B)  (A-C)Fa/CB](a*2ft)a  =  E* 

D5=  (abaTlc)[  (Ke/B)  ♦(«♦-/€]  ]  ♦  2t  *  A-B)<A-C)  pa/CB  ]  (ab*) 

♦  2(ba*2fha)  (ata)  ♦  2[ 4A-B>  (A-C) p2/CE  ]  (ba*2£ba)  (a^2f b) 

♦  (a*2fb)  (ataTJc)[  (Kr'rRv/C)*(Kr9Ke/E)  ]  =  D51  ♦  B52  pa 


D^=  (abailg)[  (Ky/C)*(Ke/E)  ]4a*2£b) 

♦  2[ (A-B)  <A-C) p2/CB](a*2£b)  lab*) 

♦  C  <Kr.vKy/C)  ♦  (Kre  Ke/E)  ]  (abailc )  (ba*2£ba) 

♦  [  (A-B)  (A-C)  pa/CE](ta»2fba)a  ♦  {ata)a  =  ♦  E^gP* 


D  3  =  [  (Ky/C)  ♦  (K9/B)  ]  (ab*TJc  )  (ba^2£ba) 

♦  2[  (A-B) (A-C) pa/CB](aba)  (ba»2£ba) 

♦  C(KryK-v/C)  ♦  (IVe  Ke/E)  ]({abaiJc  )ab2)  =  D  31  ♦  E  32pa 


=  [  (Km-/C)  ♦(K^/E)  ](aba)  (abail^) 

♦  [  (A-B)  (A-C)  pa/CE](ata)2 

♦  (Kj.^K^Kr^  Re/BC)  (Tlcaba)a  =  ♦  D22  P* 

Di  =  Re  (abaTlg/CB)  aRpy  Ky  ♦  Kr, *  e(abaTl^/CE)  ZR^ 


°0  ’  t  (abail  )  a]/CB 
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llie  denoaioatcr  of  Eg.  (3.11),  when  set  equal  to  zero, 
defines  the  chacacteristic  equation  of  the  systea. 
Exaiinatioa  of  Eq.  (3.12)  (in  wbich  Eq.  (3.11)  defines  the 
coefficients)  shews  that  as  the  roll  cate  is  increased  fron 
zero,  the  roots  of  the  characteristic  equation  will  nove  as 
a  functicn  of  the  squared  cell  rate.  Foe  analysis,  the 
region  of  operation  will  he  restricted  to  the  heginning  and 
end  of  each  booster  stage.  At  these  instants,  the  chacac¬ 
teristic  equation  is  as  given  in  Eg.  (3.12)  with  nuoerical 
values  for  the  coefficients  given  in  Table  3-1.  The  end  of 
stage  three  is  not  considered  because  the  aissile  is  no 
longer  perforaing  powered  flight  at  the  end  of  the  final 
stage. 


s»o  ♦  Cgs*  ♦  (D8i*D  a2P*)s«  ♦  (Oy^  ^2?*)  s'  (3.12) 

♦  (D6l  ♦D62P*)s‘  ♦  (D51  *052  p*)ss  ♦  (D^^  +0^2  F*)®* 

♦  (D  31^052  *  (®21*^2  P*)®*  *  D  ^S  ♦  *  0 


The  '*aigration"  of  the  systea's  roots  as  the  cell  cate 
increases  froa  zero  is  shown  in  Figs.  3-3  tfacougfa  3-12. 
Figure  3-4,  for  exaaple,  shows  that  two  coots  of  the  systea 
aove  into  the  right-half  plane  at  a  roll  rate  of  appeozi- 
aately  12.2  radiar.s/second.  Since  the  linearized  chacacter¬ 
istic  equation  was  used  to  construct  the  root  loci  shown  in 
Figs.  3-3  through  3-12,  the  precise  value  cf  roll  rate  which 
causes  instability  cannot  be  detecained  fees  the  figures. 
Shat  the  figures  point  out  is  the  general  aigeation  cf  two 
systea  coots  into  the  right-half  plane  in  all  phases  of 
operation  as  the  roll  rate  is  increased. 
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B.  BiriCT  OP  BOIl  01  SISIBB  PBBPOBBIICl 

As  shovn^  the  intcodactioD  of  roll  into  the  pitch*;aH 
controller  causes  a  ligration  of  two  sjstea  roots  into  the 
cight'>hal£  of  the  S'*plane.  As  the  roots  approach  the  iiagi- 
nary  axis  (as  the  coll  rate  is  increased)  the  systei  tegins 
to  show  sore  transient  oscillation  due  to  the  decreased 
dasping  ratio,  Ihe  settling  tise  of  the  systea  is  also 
increased.  Since  the  systea  contains  the  nonlinear  satura¬ 
tion  effect  of  the  nozzle  actuator  and  is  of  such  high  order 
(tenth),  not  such  sore  than  this  gualitatiwe  ohservation  was 
done  analytically. 

lo  study  the  effect  of  coll  on  systea  perfocaance,  a 
secies  of  siaulation  studies  was  conducted  with  the  satura¬ 
tion  effect  present.  ihe  siaulation  studies  were  peeferaed 
by  subjecting  the  systea  to  three  sets  of  step-input 
coaaands  and  constant  coll  rates,  the  first  study  consisted 
of  giving  the  systea  a  step-pitch  input  and  a  zero-yaw  input 
and  observing  the  transient  response  at  six  coll  cates 
between  2.5  and  IS  cadians/second.  Ihe  results  of  these 
siaulaticns  are  plotted  in  Figs.  3-1J  through  3-18.  The 
second  study  involved  subjecting  the  systea  to  a  zero-pitch 
input  and  a  step-yaw  input  and  conducting  the  saae  tests  as 
the  first  study.  The  results  of  the  second  study  are 
plotted  in  Figs.  3-19  through  3-24.  Finally,  the  systea  was 
given  both  a  step-pitch  coaaand  and  a  step-yaw  ccaaand  and 
the  tests  were  repeated.  Ihe  responses  of  the  systea  to  the 
conditions  of  the  third  test  are  shown  in  Figs.  3-25  through 
3-30.  Unlike  the  non-rolling  siaulations  which  were  cun  at 
separate  booster  stages  (Figs.  2-9  through  2-11),  the 
tolling  siaulaticns  were  run  over  several  consecutive 
booster  stages.  This  is  because  the  settling  tiae  cf  the 
rolling  aissile  systea  was  often  greater  than  the  turn  tiae 
of  a  single  stage. 
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Fig.  3-13  Pitch-Step  Response  at  Boll  Bate  of  2.5  cad/sec 

Figures  3-13  through  3-18  show  the  systen's  pitch  and 
yaw  responses  to  a  pitch-step  input  and  a  2ero-yaw  input  at 
the  specified  roll  cate.  Rote  the  draaatically  increased 
overshoot,  oscillation  and  settling  tine  ccnpared  to  the  the 
non-rolling  sinulations.  Figs.  2-9  thcougb  2-11.  Figure 
3-18  shows  shows  the  unstable  response  obtained  at  a  roll 
cate  equal  of  15  rad/sec. 


51 


RESPONSE  (RADIANS) 


TIME  (SEC) 


Fig.  3-16  Pitch-Step  Response  at  Boll  Bate  of  10-0  cad/sec 
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Fig.  3-19  law-step  Response  at  Boll  Bate  of  2.5  rad/sec 

Figures  3-19  through  3-24  show  the  sjstei's  pitch  and 
yaw  responses  to  a  yaw-step  input  and  a  pitch-zero  input  at 
the  specified  roll  rate.  In  these  siaulations,  as  in  the 
previous  ones,  the  overshoot  and  settling  tiie  are  increased 
above  that  observed  in  the  ncn-rolliog  siiulations  of  Figs. 
2-9  through  2-11.  Figure  3-24  shows  the  unstable  response 
at  a  roll  rate  of  15  rad/sec. 
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Fig.  J-23  law-step  Besponse  at  Boll  Bate  of  12.5  rad/sec 
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Fig.  3-25  Pitch-  8  Xaw-Step  Besponse  Boll  Bate  2.5  cad/sec 

Figures  3-25  through  3-30  show  the  systeo's  pitch  and 
yaw  responses  to  pitch-  and  yaw-step  inputs  at  the  specified 
roll  rate.  Under  these  conditions,  the  under-danped  oscil¬ 
lation  observed  in  the  two  previous  test  is  even  sore 
pronounced. 


Ihe  siBuIatioa  studies  sJiaued  that  the  systes  hegios  to 
reach  instability  at  a  roll  rate  of  approzisately  12 
radians/second.  However,  at  roll  cates  as  low  as  2.5 
radians/secoad  the  systea*s  perforaance  has  bees  severely 
degraded  froa  the  non-rolling  case  Uigs.  2-9  through  2-11). 
In  all  cases,  the  Baxiauo  peak  overshoot  is  increased  by  a 
factor  of  at  least  five  and  the  settling  tiae  by  a  factor  of 
at  least  eight.  Ihe  systea's  perforaance  is  approziaately 
that  of  an  under-daaped  oscillator.  This  response  is  caused 
by  the  two  roots  which  have  aoved  close  to  the  iaaginary 
axis  as  the  roll  rate  was  increased  and  are  doainating  the 
systea's  transient  response.  In  the  next  section,  a  coapen- 
sation  scheae  is  proposed  which  increases  the  daaping  if  the 
systea  and  partially  cancels  the  cross-coupling  acceleration 
teras  introduced  ty  roll. 


IV.  COHPgMSAlIQH  SCaj!!| 


In  the  preceding  chapter  it  was  noted  that  rolling  the 
missile  had  the  effect  of  introducing  cross'coupling  terms 
into  the  pitch  and  yaw  angular  acceleration  equations  <£gs. 
(3.6  and  (3.7)).  In  order  to  completely  restore  the 
syst6s*s  performance  to  its  non-rolling  condition,  it  is 
necessary  to  eliminate  or  cancel  out  the  cross-coupling 
terns.  Equation  (3.10)  can  fce  rearranged  into  the  form: 
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The  left  side  of  £q.  (4. 1)  contains  the  "dynamics"  of  the 
system.  In  other  vords,  it  gives  the  acceleration  and 
velocity  terns  which  determine  the  angular  acceleration  in 
the  pitch  and  yaw  channels.  The  only  way  to  eliminate  the 
cross-coupling  terms  directly  would  be  to  make  changes  to 
the  left  side  of  Eg.  (4.1)  which  would  cancel  the  terns,  Qs 
and  Bs6.  This  would  involve  dramatic  changes  to  the  missile 
body,  thrust  vector  control  system  and  flight  performance 
and,  therefore,  will  not  be  attempted. 

The  right  side  of  Eg.  (4.1)  constitutes  the  "elec¬ 
tronics"  of  the  system.  It  contains  a  sun  of  those  terns 
which  are  added  at  the  feedback  junction  then  passed  through 
the  nozsle  actuator.  This  side  of  the  equation  can  be 
changed  easily  by  altering  the  electronic  network  that  makes 
up  the  feedback  loop. 


A.  CBOSS-CHAIMBl  FBBOBACB 

Assuae  that  the  pitch  cate  and  yav  cate  cao  be  acca- 
cately  aeasaced  oc  estiaated  at  any  tiae  daciog  flight. 
Also«  assuae  that  the  tecas,  G;  fl;  C  and  B,  all  functions  of 
systea  pacaaetecs,  can  be  aeasuced  oc  estiaated  accucately. 
Then,  addition  of  the  feedback  tecas  of  pitch-cate-to-yaw- 
input  and  yaa-rate-to-pitch-input  pecait  cancellation  of  the 
ccoss-coupling  tecas.  The  block  diagcaa  cealization  of  the 
pcoposed  change  is  shown  in  Fig.  4-1.  Equation  (4.2)  gives 
the  coapensated  systea*s  tcansfocaed  diffecential  equation 
in  aatcix  foca. 
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Equation  (4.2)  can  be  ceduced  to  Eg  (4.3)  by  canceling 
the  ccoss-coupling  tecas,  Cs  and  BsO,  fcoa  both  sides  of 
the  equation.  Coipacison  of  Eq.  (4.3)  with  Eq.  (3.5)  shows 
that  the  tcansfocaed  diffecential  equation  of  the  colling 
systea  has  been  ceduced  to  its  ociginal  non-colling  foca. 
Note  that  all  that  was  done  in  this  coapensation  scheae  was 
to  subtcact  the  ccoss-coupling  teen  fcoa  the  coanand  input 
to  the  nozzle  actuatoc  so  than  when  it  is  added  back  on 
downstceaa  the  cesult  will  be  no  effective  ccoss-channel 
teens  into  eithec  ccntcol  path. 


Rre  S0  -  0) 
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(4.3) 


UDfoctunately ,  the  pcoblee  is  not  as  sisple  as  £q.  (4.2) 
iaplies.  The  laio  ptobles  is  that  the  aaount  of  control 
nhich  can  be  used  to  force  the  systea's  response  is  United 
by  the  saxiaan  deflection  angle  of  the  exhaust  nozsle.  This 
nonlinear  saturation  of  the  actuator  is  not  present  in  Eg. 
(4.2).  Figure  4-2  suanarizes  the  problea. 


Fig.  4-2  Plant  Dynaaics  and  Control  Input  liaitations 

Figure  4-2  shows  that  no  natter  what  feedback  scheae  is 
enployed  to  generate  the  pitch  and  yaw  angular  acceleration 
control  signals,  they  ace  United  and,  therefore,  so  is  the 


opecation  over  which  the  pecfociance  of  the  sjstea 
can  be  cootcolled  and  the  degree  of  control  possible. 

Figure  4-3  shoes  the  cross-coupling  teres  and  the 
atteepted  cancellation  of  thee  by  the  cross-channel  feedback 
teres  under  the  coepensation  scheee  just  proposed.  Ihe 
sieulation  shown  was  perfoceed  during  the  first  stage  of 
operation  with  the  systee  subject  to  a  pitch-step  coeeand 
and  a  eero-yae  coeeand.  Note  bow  the  eagnitude  of  the 
cross-channel  feedback  tern  is  effectively  saturated  by  the 
nozzle  actuator  and  cannot  follow  the  cross-coupling  tern 
outside  a  seall  region  of  operation.  Ihe  sieulation  was 
made  by  subtracting  the  error  tern  generated  by  the  differ¬ 
ence  in  coeeand  and  feedback  so  that  only  the  cross-channel 
feedback  tern  is  fed-through  the  nozzle  actuator. 

B.  NOMLIHBAB  FEEOEICR 

Exaeination  of  Figs.  2-5,  2—9,  2-10  and  2-11  shows  that 
the  tine  the  nozzle  actuator  takes  to  reach  eaxieue  deflec¬ 
tion  when  given  a  large  input  coeeand  (>>  eaxieue  deflec¬ 
tion)  is  insignificant  (0.05  sec)  when  ccnpared  to  the 
settling  tine  of  the  entire  pitch-yaw  control  systee.  Thus, 
the  actuator  is  functioning  aleost  like  a  switch  in  all 
regions  of  operation  except  near  the  coeeand  input  (steady 
state).  For  exaeple.  Fig.  4-4  shows  the  nozzle  deflection 
in  the  pitch  plane  during  a  sieulation  when  the  systee  was 
subject  to  a  pitch-step  input  and  a  roll  rate  of  2.5 
radians/second.  The  conclusicn  to  be  drawn  fron  Fig.  4-4  is 
that  since  the  aeplitude  of  the  effective  control  input  to 
the  plant  cannot  be  increased  above  saturation,  the  only 
alternative  is  to  "throw  the  switch'*  a  little  sooner. 

The  cross-channel  feedback  coepensation  scheee  just 
proposed  will  cause  a  reduction  in  the  effective  servo 
coeeand  and  result  in  earlier  reeoval  of  the  eaxieue  nozzle 
deflection  (turning  off  the  switch).  However,  this  approach 
can  at  best  only  fallow  the  cross-coupling  terns  and  can  not 
"anticipate"  the  need  to  switch  the  deflection  angle  froe 


positive  naxiaun  to  negative  tiaxiaua.  Another  aethod  to 
acconplisb  the  objective  o£  svitching  the  nozzle  earlier  is 
to  var;  the  velocity  feedback  gains,  and  . 
Specifically,  it  is  desired  that  the  velocity  feedback  gains 
be  low  daring  the  initial  phases  of  operation  when  the  error 
between  coaaand  and  position  is  large  and  increase  as  the 
systen  approaches  the  coaaandcd  value.  Ihis  will  result  in 
the  sane  basic  systen  operation  that  was  present  tefore 
conpensation  (doainant  coaplex*pair  oscillation)  but  with 
low  danping  when  the  error  is  high  and  higher  danping  as  the 
systen  approaches  steady  state. 

The  nonlinearity  of  the  system  prevents  the  use  of  nost 
analytical  tools  which  night  be  used  to  deteraine  algorithns 
for  R  pe  *  Instead,  a  trial  and  error  approach  was 
used  based  on  the  observations  of  the  systea*s  perfornance 
during  earlier  siaulations  and  siaulation  studies  perforned 
specifically  to  deteraine  appropriate  paraaetecs. 

Assune  that  the  algorithn  for  deternining  the  velocity 
feedback  gain  consists  of  a  constant  which  is  reduced  by  an 
amount  proportional  to  the  error  between  coaaand  and  posi¬ 
tion.  Ihis  will  result  in  low  danping  if  the  error  is  high 
and  increased  danping  as  the  error  decreases.  For  exaaple, 
if  V  and  H  are  ccnstants,  the  pitch-channel  velocity  feed¬ 
back  gain  can  be  determined  from: 


=  V  -  «  0^  -  0 


(4.4) 


Furthernore,  the  siaulations  shown  in  Figs.  3-13  through 
3-30  indicated  that  the  systen  becoaes  nore  undeidanped 
(roots  migrate  toward  the  iaaginary  aiis)  as  the  roll  rate 
is  increased.  So,  the  velocity  feedback  gains  were  nade 
proportional  to  the  roll  rate.  Nornalizing  the  algorithn  so 
that  the  relationship  does  not  becone  zero  if  the  roll  rate 
is  zero  yields  Egs.  (4.5)  and  (4.6). 


Pignre  4-5  shows  sanple  cesults  of  sisulation  studies 
which  were  conducted  foe  various  values  of  V  and  R  and  over 
various  stages  of  operation  and  roll  rates.  Ihe  sinulations 
showed  that  acceptable  perforaance  was  obtained  for  V  s  4 
and  R  -  3. 

Figures  4-6  through  4-14  show  simulaticn  results  which 
were  nade  b;  using  the  cross-channel  feedback  and  variable 
velocity  feedback  gain  conpensation  schenes  together.  For 
coaparison*  the  uncoapensated  responses  are  also  shown.  Ihe 
coapensated  siaulation  responses  were  obtained  under  the 
sane  operating  conditions  as  those  shown  in  Figs.  3-13 
through  3-30. 
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Fig.  4-5  Saaple  Algoritha  Eesponse  lest  Siaulation 


RESPONSE  (RADIANS) 


NSE  (RADIANS) 


Uncompensated  Pitch 
Compensated  Pitch 
Uncompensated  Yaw 


Compensated  Yaw 


I 


^Uncompensated  Yaw 
/Compensated  Yaw 
/Uncompensated  Pitch 
Compensated  Pitch 


RESPONSE  (RADIANS) 


RESPONSE  (RADIANS) 


RESPONSE  (RADIANS) 


RESPONSE  (RADIANS) 


TIME  (SEO 


Fig  4-14  Pitch-  fi  Yaw-Step  Besponse  Bell  Bate  7.5  cad/sec 


The  siaulatioDS  indicate  that  substantial  pecforaance  resto¬ 
ration  is  possible  with  the  proposed  coipensation  scheies. 
However,  the  restoration  is  not  conplete  and  is  still  luch 
worse,  in  terns  c£  overshoot  and  settling  tine,  than  the 
non-rolling  nissile. 


V.  COICIOSIOHS 


The  coBpensation  scheae  proposed  partially  restores  the 
pre-roll  pecforaaDce  of  the  aissile.  Table  5-1  suaaarizes 
the  level  of  perfcraance  degradation  caused  by  roll  and  the 
restoration  ehich  was  obtained  vith  this  ccapensation  schene 
under  pitch-  and  yaa-step  inputs  and  a  roll  rate  of  2.5 
rad/sec.  As  iadicated,  the  degradation  is  significant  even 
at  this  relatively  low  (and  stable)  roll  rate.  The  coapen- 
sated  systea's  perforaance  is  still  far  votse  than  that  of 
the  non-rolling  systea. 

The  coapensation  scheae  relies  priaarily  on  the  nonli¬ 
near  daaping  effect  of  the  variable  velocity  feedback  gains 
which  "anticipate"  the  approach  of  the  coaaanded  position. 
Since  the  cross-coupling  terns  introduced  by  roll  are  not 
directly  cancelled  by  the  proposed  aodif ication,  the  region 
of  stable  operation  and  the  degree  to  which  the  systea*s 
perforaance  can  be  restored  is  directly  controlled  by  the 
roll  rate.  Purtheraore,  the  saturation  caused  by  the  nozzle 
deflection  will  prevent  any  coapensation  scheae,  which  only 
aodifies  the  control  input,  froa  fully  restoring  the 
systeB*s  perforaance.  To  coapletely  restore  the  systea  to 
its  pre-roll  perforaance  at  any  roll  rate,  additional 
control  inputs  to  the  systea  plant  dynaaics  (such  as,  side 
thrusters  or  additional  vectoring  of  the  main  rocket  aotor) 
are  necessary. 
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The  DSI/360  prograo  used  to  sinulate  the  Doo'colling 
■issile  is  shovD  beloM.  The  pcograa  consists  of  three 
stages*  initial*  derivative  and  teininal.  The  initial 
segnent  is  used  to  define  constants  and  calculate  functions 
of  the  constants  uhich  will  net  change  throughout  the  sinu- 
laticn  run*  The  derivative  section  is  where  the  actual 
integrations  reguired  to  perforo  the  siaolations  are 
conducted.  This  progran  is  constructed  to  follow  the  block 
diagraw  of  Fig.  2-1  directly.  Mote  that  the  nozzle  actuator 
transfer  function  was  broken  down  into  a  three-step  integra¬ 
tion.  The  terninal  stage  determines  the  end  condition  of 
the  sinulation  run  and  the  output  data.  This  data  was 
formed  into  a  data  file  and  a  packaged  graphics  program 
(DISSELA)  was  used  to  construct  the  figures  shown  in  the 
teit.  The  parameter  values  listed  in  this  example  program 
are  those  during  stage  one. 


*  miTIAlIZE  COBSTAIIIS 

CONST  TST=7496000*L=J63.J1=11048.8. KTH=1.0,KBTH=1.75*. 
^  «A=150*WB=15,PB=0.9*FATE=9  30.0*m=  189140 

*  COBPDTE  SET  TALOES 
INITIAL 
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?=INTGBL  (0,E1D01) 


«  OOTPOI  AND  CONIfiOl  STAIE8EN1S 
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SAVE  0.015,P 
PBINT  0.015,P 
END 
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Unlike  the  non-rolling  systea  sinulations,  the  rolling 
■issile  sinulations  were  required  to  run  over  several 
consecutive  booster  stages,  this  required  the  addition  of  a 
'dynaniC  segaent  of  siaulatioo.  Ihe  *1F*  stateaents,  shown 
in  the  exaaple  siaulatioo  below,  are  used  to  switch  froa  one 
set  of  paraaeters  to  another  as  one  stage  burns  off  and 
another  starts.  The  switching  is  controlled  by  the  elapsed 
tine  of  the  siaulatioo  run  coapared  to  the  burn  tiae  of  each 
stage. 
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KTH=1.0,KeTH=1.  5,ilA=15O,«B=15,PB=0.S,RClt=7.5 
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INITIAL 
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EESET=0.0 

^INITIALIZE  FIBST 
TST=7496000 
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J1=110«8.8 
02=136,  1 
EAIE=930,0 
H1=189140 

* 

*  DECIDE  WHICH  STAGE 

DINAHIC 

*  ENTER  STAGE  TWO 

IF  (TIHE.LT. 150) 
TST=1125000 
I  =  230 
01=4476.4 
BA1E=81 .49 
nl=40340 
BBSEI=150 

*  ENTEB  STAGE  THBEE 

IF  (TIHE.LT. 509) 
TST=224900 
1=100 
01=858.3 
02=50.0 
EATE=14.75 
nl=8140 
BESET=509 

10  CONTINUE 
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*  EBBOB  SIGM&IS  FOSlllON  FEEDBACK 

£P  *  PCOM-F 
£1  =  ICOH-Y 

*  EBBOB  SIGNALS  VELOCITY  FEEDBACK 

£IP  =  EP-IKBTH*P1001) 

ElY  =  EY-iKBTfl*Y1001) 

CMP  =  KTH*EIP 
CBY  -  RTB«EIY 

*  SEBVO  TBAMSFEB  FUNCTION  IN  FITCH  PATH 

SVF  =  IMTGBLIO.SVPDOT) 

SVPOOT  3  INT6BI  lO.SVF^DT) 

SVF2DT  *  IHTGBL  |0-SVP3DTi 

SVP3DT  =  A*(CMP-S7P)  -  E*SVP2DT  -  C*SiPDOT 
«  PITCH  DEFLECTION 

POEFL  »  LIMIT (-0.0873, 0. 0873, S7P) 

*  SEBVO  TBAMSFEB  EUMCTIOM  IN  YAM  PATH 


SVY  =  IMTGBL (O.SVYDOT) 


SVYDOT  =  INTGBL  (0,SVY20T} 
SVI2DT  -  IMTGBL  {0,SVY3DT) 
SVY3DT  =  A*ICMI-SVY1  -  E* 
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SVI3DT  =  A*(CMI-SVY)  -  B*SVY2DT  - 

*  YAH  DEFLECTION 

YDEFL  =  LIMIT  (-0.0873, 0.0873, SVY) 

JY  =  J1*aASS 

JZ  »  JY 

JX  -  J2«BASS 

ZP  3  TST*L/JY 

2Y  *  TS1*L/JZ 

*  CBOSS  COUPLE  MULTIPLISBS 

ALPHA  =  IJZ-JX) ♦BOLL/JY 
BETA  -  iJl-JY)*BOLL/JZ 
«  CBOSS  COUPLE  TEBH  YAH  TO  PITCH 
CCYP  s  ALPHA*I100T 

*  CBOSS  COUPLE  TEEM  PITCH  TO  YAH 

CCPY  *  BBTAtPlDOT 

*  PITCH  ANGULAB  ACCELEBATIOM 

P200T  (ZP*FDEFLI  *  CCYP 

*  YAH  ANGULAB  ACCELEBATlOM 

Y2DOT  =  (ZY*YDEFL)  ♦  CCPY 

*  PITCH  ANGULAB  VELOCITY 

P100T  s  INIGBL jO,P2DOI) 

P  =  INIGBL  10, PiDOT) 

*  YAH  ANGULAB  VELCCITY 

Y100T  -  IMTGEL10,Y200T) 

^  Y  »  INTGBL  (0,11DOT) 

*  OUTPOT  AND  COMIBOL 
TEBMINAI 
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